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ABSTRACT 

We are undertaking a program to measure the characteristics of the intracluster light (total flux, 
profile, color, and substructure) in a sample of 10 galaxy clusters with a range of cluster mass, 
morphology, and redshift. We present here the methods and results for the first cluster in that 
sample, A3888. We have identified an intracluster light (ICL) component in A3888 in V and r that 
contains 13 ± 5% of the total cluster light and extends to TOOh^kpc (~ 0.3r2oo) from the center of 
the cluster. The ICL color in our smallest radial bin is V — r = 0.3 ± 0.1, similar to the central cluster 
ellipticals. The ICL is redder than the galaxies at 400 < r < TOOhf^kpc although the uncertainty in 
any one radial bin is high. Based on a comparison of V — r color with simple stellar models, the ICL 
contains a component which formed more than 7 Gyr ago (at z > 1), coupled with a high metallicity 
(I.OZq < Zjcl ^ 2.5Zq), and a more centralized component which contains stars formed within the 
past 5 Gyr (at z ~ 1). The profile of the ICL can be roughly fit by a shallow exponential in the outer 
regions and a steeper exponential in the central region. We also find a concentration of diffuse light 
around a small group of galaxies l^h^Mpc from the center of the cluster. In addition, we find 3 low 
surface brightness features near the cluster center which are blue (V — r = 0.0) and contain a total 
flux of 0.1M*. Based on these observations and X-ray and galaxy morphology, we suggest that this 
cluster is entering a phase of significant merging of galaxy groups in the core, whereupon we expect 
the ICL fraction to grow significantly with the formation of a cD galaxy as well as the in-fall of groups. 
Subject headings: galaxies: clusters: individual (A3888) — galaxies: evolution — galaxies: interactions 
— galaxies: photometry — cosmology: observations 



1. INTRODUCTION 

Galaxy clusters contain a population of stars which 
are not members of individual galaxies but which are 
bound to the cluster potential, producing diffuse in- 
tracluster light (ICL). This ICL component has been 
found in a number of clusters through surface bright- 
ness measurements and direct detections of resolved 
stars including planetary nebulae, red giants, super- 
novae, novae, and globular cluster systems. These in- 
vestigations indicate that the optical ICL comprises be- 
tween 5% and 50% of the total optical cluster lumi- 
nosity (see IFeldmeier et all 120041 IGonzalez et~all 120051: 
IZibetti et alJl2005l and references therein). Conclusions 
on the color of the ICL vary wi dely from blue to red, with 
and without color gradients llSc hombertl Il988t [Mackie 
ITflfll IGonzalez et ail fZcM IZibetti et a,l .1 120051) . Cur- 
rent measurements of the shape of the ICL generally 
favor a double deVaucouleurs profile such that one func- 
tion fits the BC G and the second function fits the ex- 
tended envelope (IGonzalez et al.ll2000t iBernstein et alJ 
119951 IZibetti et all 120051) . Examples of tidal features 
such as plumes and bridges are found in multiple clus- 
ters as evidence of interactions which add stars to the 
ICL lIGregg fc Wesdll998t ICalcaneo-Roldan et alJl200?j 
IFeldmeier et al.M2004|) . Long-slit spectroscopy of A2199 
shows that the intracluster stars there h ave the same 
velocity dispersion as the cluster galaxies ([Kelson et al.l 
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2002), confirming that, in at least one cluster, the intr- 
acluster stars are not bound to the individual galaxies 
but trace the overall cluster potential. Conversely, intr- 
acluster planetary nebulae studies show evidence for less 
relax ed velocities ([Arnaboldi et al J [20041: iGerhard et alJ 
2005). There is no consensus on the velocity distribution 
of intracluster stars. 

The intracluster light is a fossil remnant of cluster for- 
mation and evolution and can be used to study the dom- 
inant physical processes involved in galaxy evolution in 
clusters. Hierarchical dark matter simulations suggest 
that galaxies falling into dense regions would loose most 
of their mass. When mechanisms such as radiative cool- 
ing and star formation are included in the simulations, 
galaxies which are composed of a central dense core of 
stars do retain most of their stellar mass throughout clus- 
ter in-fall, but lose some stars to the cluster potential. 
State of the art simulations are able to predict the ex- 
istence of this intracluster population, but basic ques- 
tions as to its properties can only be answered by un- 
derstanding which physical mechanisms are important. 
This work seeks to answer when and how intracluster 
stars are formed by studying the total flux, profile shape, 
color, and substructure in the ICL as a function of cluster 
mass, morphology, and redshift. 

Observations of the total flux in the ICL over a sample 
of clusters will allow us to identify the effects of clus- 
ter environment on galaxy evolution. For example, a 
high mass cluster should have a higher ICL fraction than 
low-mass clusters if ram pressure stripping or harassment 
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are dominan t mech anis ms. In fact, simulation s by both 
iLin fc Mohrl ((1)04) and lMurante et alJ l|2004D predict a 
strong relation between ICL fraction and mass. If, how- 
ever, galaxy-galaxy merging is the dominant mechanism, 
and most of the galaxy evolution happens early on in 
cluster collapse, then the ICL should not correlate di- 
rectly with current cluster mass. The existence of a cD 
galaxy in a cluster is evidence of a rich merger history 
and therefore morphology should also correlate strongly 
with ICL fraction. The ICL fraction will be affected by 
rcdshift, since with time comes an increased number of 
interactions. 

Observations of the color and substructure of the ICL 
will help to identify the origin, formation epoch, metal- 
licity, and possibly progenitor morphologies of cluster 
galaxies. For example, if the ICL is as red or redder than 
the bright cluster ellipticals, it is likely to be a remnant 
from the early epochs of cluster formation with little re- 
cent accretion of tidally disrupted systems. If the ICL is 
bluer than the galaxies, then some recent accretion has 
occurred, either from ellipticals with low metallicity or 
spirals with younger stellar populations. While multiple 
mechanisms are likely to play a role in the complicated 
process of formation and evolution of clusters, important 
constraints can come from ICL measurement in clusters 
with a wide range of properties. 

In addition to constraining galaxy evolution, the ICL 
is an important baryonic component in clusters. The 
ICL, which is typically not included in the baryon cen- 
sus, will contribute to the baryon fraction of clusters and 
that contribution is likely to change with time. If the 
ICL fraction does indeed evolve with redshift and is a 
significant fraction of the total cluster light, it will sys- 
tematically bias the inferred red shift depend e nce o f the 
baryon fraction. Recent work bv lAllen et alJ l)2004(l has 
used a change in baryon fraction with redshift of only a 
few percent to constrain cosmological parameters. When 
doing such precision cosmology it will be necessary to in- 
clude ICL in the cluster light budget. 

The ICL may also play an important role in the 
global properties of the intracluster medium (ICM). It 
has recently been suggested that an intracluster stel- 
lar population (ICSP) can account for at least some 
amount of heating an d metal enrichment of the intra- 
cluster medium (ICM) jZaritskv et alJ2004tlLin k. Mohrl 
120041 iDomainko et al.ll2004[) . Considering only super- 
novae (SNe) within galaxies, the ful l metallicity of th e 
ICM can not easily be accounted for ijLin fc Mohrl l2004). 
However, since intracluster supernovae are in situ in the 
ICM, they contribute directly to the metallicity of the 
ICM, and will have therefore a direct impact on its abun- 
dance. Although these authors find that the ICL can- 
not fully account for the high abundance of the ICM 
(~ 0.3Zq), further studies are warranted to quantify just 
how many ICSNe there are. Even if the ICSP can not 
account for the full metallicity of the ICM, it is possi- 
ble that this population is responsible for the metallicity 
gradient found in clusters. If true, a correlation between 
ICL flux and abundance gradients in clusters should ex- 
ist. 

In this paper we present the methods of this survey as 
well as measurements of the color, total flux, and profile 
shape for the first cluster in our sample, A3888. In 32] 
we discuss characteristics of the entire sample. Details of 



the observations and reduction are presented in and 
21 including flat-fielding and sky background subtraction 
methods. In §21 we discuss object detection and removal 
as well as cluster membership. In §15] we describe results 
followed by a discussion of accuracy limits in §3 In JHl 
we present a discussion of the results. Conclusions are 
summarized in §5] 

Throughout this paper we use H = 70km/s/Mpc, CIm 
— 0.3, flA — 0.7 which gives 3.5 kpc/arcsecond at the 
distance of A3888. 

2. THE SAMPLE 

The ten galaxy clusters in our survey were selected to 
meet two general criteria. First, each cluster has a pub- 
lished X-ray luminosity which guarantees the presence of 
a cluster and provides an estimate of the cluster's mass. 
Second, all are at high galactic and ecliptic latitude along 
lines of sight with low HI column density. This minimizes 
complications due to scattered light from galactic stars 
and zodiacal dust and from variable extinction across the 
cluster field. Of the clusters that meet the above qual- 
ifications, we selected ten clusters as the beginning of a 
statistical sample which is representative of a wide range 
in cluster characteristics, namely redshift, morphology, 
spatial projected density (richness), and X-ray luminos- 
ity (mass). 

To the extent possible, we also selected clusters for 
which mass estimates and membership information are 
available in the literature. For example, in addition to 
published X-ray masses, three of the clusters have mass 
estimates from gravitational lensing measurements. Pub- 
lished redshift surveys provide velocity dispersions and 
membership information for all but 2 clusters in the sur- 
vey. Those 2 clusters, as well as 5 others with small 
numbers of published velocities, were included in a red- 
shift survey we undertook with IMACS on Magellan I 
(Baade) . With these additional observations the physical 
properties of all clusters in our sample can be compared 
to the ICL characteristics. Table 1 lists the relevant in- 
formation for the cluster sample. 

The sample is divided into a "low" (0.05 < z < 0.1) 
and "high" (0.15 < z < 0.3) rcdshift range which we 
have observed with the 1 meter Swope and 2.5 meter 
Du Pont telescope respectively. The bottom end of the 
redshift range is limited by the field of view of the 1 
meter telescope and detector , which corresponds to 0.9 x 
1.4/ifo 1 Mpc (14.8' x 22.8') at z = 0.05. This field of 
view allows us to measure the ICL as well as off-cluster 
background flux in the same image for all clusters in 
the sample. The top end of the redshift range reflects 
X-ray data availability and the increasing difficulty of 
measuring diffuse sources at high redshift due to (1 + z) A 
surface brightness dimming. 

2.1. A3888 

This paper focuses on one cluster in our sample, 
A388 8, which is a richness class 2 cluster at z=0.15 
(jAbell et all 11989(1 . This Bautz- Morgan type I-II clus- 
ter has no cD galaxy; instead the core is comprised of 3 
distinct sub-clumps of multiple galaxies each. At least 
2 gal axies in each of these clumps are confir med mem- 
bers <|Teaeue et al.lll990tlFimbblet et aT]|2002j). On large 
scales (286, 535, and 714 /i^kpc) iGirardi et all <|1997f) 
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find a unimodal distribution for this cluster with no de- 
tected substructures in either the galaxy spatial or ve- 
locity distribution. The projected spatial distribution of 
galaxies in A3888 is slightly elon gated with an ellipticity 
of 0.43 l|Struble fc Ftaclasill994l) . X-ray surface bright- 
ness from XMM observations also indicate an elongated, 
single-peaked distribution for the hot gas. The cluster 
contains an X-ray bright Seyfert I galaxy located at a 
projected distance of roughly 6 00h^ n kpc from the clus- 
ter center ijReimers et al.lll996j) . 

The mass of A3888 can be estimated from two dif- 
ferent sets of observations. iReinrich fc Bohringerl ( 2002) 
calculate gravitational mass based on pointed ROSAT 
PSPC count rates and the ROSAT-ASCA Lx-Tx relation 
ijMarkevitchl 1998]) . Assuming an isothermal distribution 
and employing hydrostatic equilibrium, they find M200 
= 25.5±^°4 8 x 10 14 h-J M , where r 200 = 2.8/i 70 1 Mpc 
which is defined as the radius within which the mean 
mass density is equal to two hundred times the criti- 
cal density. In a complementary method, mass can be 
determined from published galaxy velocity dispersions. 
Based on redshifts for 50 m ember galaxies locat ed within 
a radius of S.llh^Q Mpc llTeaeue et al.lll990D and us - 
ing the method described bv lGirardi fc Mezzettil ll200l . 
we find that the mass of A3888 within r2oo is -M200 = 
40.2±3,° 4 6 x I0 u hjj Mq. For the purpose of this work, 
these two mass estimates are in good agreement, par- 
ticularly since this cluster is elongated and likely not in 
dynamic equilibrium. 

3. OBSERVATIONS 

Observations for the entire sample of 10 clusters have 
been completed. The "high" redshift observations were 
made with the du Pont 2.5m telescope at Las Campanas 
Observatory. We used the thinned, 2048 x 2048 Tek- 
tronix "Tek#5" CCD with a 3e~ /count gain and 7e~ 
readnoise. The pixel scale is 0.259"/pixel pixels), 
so that the full field of view is 8.8' on a side, correspond- 
ing to l.S/iy^Mpc per frame. Data was taken in two fil- 
ters, gunn-r (A = 6550 A) and V (A = 5400 A). These 
filters were selected to provide some color constraint on 
the stellar populations in the ICL, while avoiding flat- 
fielding difficulties at longer wavelengths and prohibitive 
sky brightness at shorter wavelengths. 

Observing runs occurred on August 19-25, 1998, 
September 2-10, 1999, and September 22-27, 2000. 
Specifically, A3888 was observed on the nights of Septem- 
ber 2 and 8, 1999 and September 22-25, 2000. Both ob- 
serving runs took place in the days leading up to new 
moon. September 2, 1999 was the only non-photometric 
night, and only 3 cluster images were taken on that night. 
These were individually tied to the photometric data. 
The average seeing during the 1999 and 2000 run was 
1.77 and 0.93 arcseconds respectively. Across both runs, 
we observed A3888 for an average of 5 hours in each 
band. In addition to the cluster frames, night sky flats 
were obtained in nearby, off-cluster, "blank" regions of 
the sky with total exposure times roughly equal to one 
third of the integration times on cluster targets. Night 
sky flats were taken in all moon conditions. Typical V— 
and r— band sky levels during the run were 21.3 and 21.1 
mag arcsec -2 , respectively. 

Cluster images were dithered by one third of the field 
of view between exposures. The large overlap from the 



dithering pattern gives us ample area for linking back- 
ground values from the neighboring cluster images. Ob- 
serving the cluster in multiple positions on the chip 
is beneficial because upon combination large-scale flat- 
fielding fluctuations will be reduced. Integration times 
were typically 900 seconds as a compromise between 
signal-to-noise and moderating the number of saturated 
stars. 

Observations of the "low" redshift clusters will be dis- 
cussed in a future paper. 

4. REDUCTION 

In order to create a single, mosaiced image of the 
cluster with a uniform background level and accurate 
resolved-source fluxes, the images must be bias and dark 
subtracted, flat-fielded, flux calibrated, background- 
subtracted, extinction corrected, and registered before 
combining. These issues are dealt with as described be- 
low. 

4.1. Bias and Dark Subtraction 

Pre-processing of the data, including overscan, bias, 
and dark subtraction, was done in the standard man- 
ner using mainly IRAF tasks. The average bias level 
was stable at ~ 800 counts, changing by 1% through- 
out the night. There is structure in the bias in the form 
of random fluctuations, as well as a highly-repeatable, 
large-scale ramping in the first 500 pixels of every row. 
To remove this structure, we first fit an 8th order poly- 
nomial to 140 overscan columns and subtract that fit, 
column by column, from each image. We further av- 
erage together ten bias frames per night with 3cr cos- 
mic ray rejection and then boxcar smooth in the vertical 
direction before subtracting from the data. We choose 
to smooth in the vertical direction because we have al- 
ready removed vertical structure in the previous process- 
ing step. Test-reduction of the bias frames themselves 
with this procedure reveals no remaining visible structure 
and each frame has a mean level of counts to within 
±0.05 counts. 

Twenty-five dark exposures were taken per observing 
run. We averaged these together with 3<r rejection to 
look for structure or significant count levels in the dark 
current. The mean dark count is 0.6 counts/900s, which 
is less than 0.08% of the sky level, and is therefore not 
significant. However, even at this small count level, there 
is some vertical structure in the dark which amounts to 1 
count/900s over the whole image. To remove this large- 
scale structure from the data images, the combined dark 
frame was median smoothed over 9x9 pixels (2.3"), 
scaled by the exposure time, and subtracted from the 
program frames. Small scale variations were not present 
in the dark. Errors in both the bias and dark subtraction 
due to structure in the residuals are an additive offset to 
the background level. These are included in our final 
error budget based on an empirical measurement of the 
stability of the background level in the final combined 
image (see JJ. 

4.2. Flat Fielding 

The accuracy of any low surface brightness measure- 
ment is limited by fluctuations on the background level. 
A major contributor to those fluctuations is the the large- 
scale flat-fielding accuracy. Pixel-to-pixel sensitivity 
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variations were corrected in all cluster and night-sky flat 
images using nightly, high S/N, median-combined dome 
flats with 70,000 - 90,000 total counts. After this step, 
a large-scale illumination pattern of order 1% remains 
across the chip. This was removed using combined night- 
sky flats of "blank" regions of the sky. To make these 
night-sky flats, objects in the individual blank sky frames 
wer e first masked before com bination. We used SExtrac- 
tor ijBertin &: Arnoutsl H996) to identify all sources with 
a minimum of 6 pixels and a total flux of 2a above the 
sky background. Mask sizes were increased by 4-7 pixels 
over the semi-major and semi-minor axes from the ob- 
ject catalogs to insure object rejection. The masked im- 
ages were then median combined with 2a rejection. This 
produced an image with no evident residual flux from 
sources and kept the large scale illumination pattern in- 
tact. Fluctuations are less than 0.1% peak to peak on 10" 
scales. The final combined night-sky flats were then me- 
dian smoothed 7x7 pixels (2"), normalized, and divided 
into the program images. The illumination pattern was 
stable among images taken during the same moon phase. 
Program images were corrected only with night sky flats 
taken in conditions of similar moon. The contribution of 
flat fielding to our final error budget is measured empir- 
ically, as described in |J3 

4.3. Non-linearity 

Although the ICL measurement is based on a low num- 
ber of counts, photometric calibration is based on bright 
standard stars. Accurate calibration is then dependent 
upon the CCD having a linear response to flux. To as- 
certain if Tek#5 was linear with flux over a wide dy- 
namic range, a consecutive chain of dome-flat images 
were taken, with exposure times of 2 - 100 seconds, cor- 
responding to approximately 300 — 15, 000 counts per 
pixel. Multiple passes through the exposure time se- 
quence (increasing and decreasing) were made to rule 
out any effects from fluctuating lamp flux. We find that 
the Tek=#=5 CCD does have an approximately 2% non- 
linearity, which we fit with a second order polynomial 
and corrected for in all the data. The same functional fit 
was found for both the 1999 and 2000 data. Note that 
exposure times used for all observations are long enough 
that shutter performance is not a problem. The uncer- 
tainty in the linearity correction is incorporated in the 
total photometric uncertainty discussed below. 

4.4. Photometric Calibration 

Photometric calibration was performe d in the usual 
manner. Fifty t o seventy standard stars j|Landoltlll992t 
I.Torgensenlll994|) were observed per night per filter over 
a range of airmasses. Stellar magnitudes were measured 
with an aperture size of 5 x the full width at half maxi- 
mum (fwhm), where the fwhm of the stars in the images 
was determined using SExtractor. We choose this aper- 
ture size as a compromise between aperture correction 
and added background noise. Photometric nights were 
analyzed together; solutions were found in each filter for 
a zero-point and extinction coefficient with an RMS of 
0.03 magnitudes (r) in Sept. 1999 and 0.02 magnitudes 
(r and V) in Sept. 2000. These uncertainties are a small 
contribution to our final error budget, but are included 
for completeness as discussed in Observing the same 
cluster field for long periods throughout the night allows 



us to measure an extinction coefficient from stars in the 
cluster fields, which we find is fully consistent with the 
extinction coefficient measured from the standard stars. 
The three exposures taken in non-photometric condi- 
tions were individually tied to the photometric data using 
roughly 10 stars well distributed around each frame to 
find the effective extinction for that frame. We find a 
standard deviation of 0.03 within each frame, with no 
spatial gradient in the residuals. 

We have compared our V— and r— band magnitudes 
for hundreds of galaxies in the cluster with R— band mag- 
nitudes from t he Las Campanas/AA T Rich Cluster Sur- 
vey (LARCS, Pimbblc t"etaD l20f)l . To the detection 
limit of the LARCS photometry, and adopting a single 
average galaxy color to convert between filters, the two 
samples are consistent with an RMS scatter of 0.07 mag- 
nitudes. 

4.5. Sky Background Subtraction 

An important issue for accurate surface brightness 
measurement is accurate identification of the background 
sky level. The off-cluster background level in any im- 
age is a combination of atmospheric emission (airglow) 
and light from extra-terrestrial sources (zodiacal light, 
moonlight, starlight, starlight scattered off of galactic 
dust). Zodiacal light comes from solar photons scattered 
off of ecliptic dust and is therefore concentrated in the 
ecliptic plane, which, along with the galactic plane we 
were careful to avoid in sample selection, so the extra- 
terrestrial background light will not vary spatially. Light 
from the extra-terrestrial sources will additionally be 
scattered into the field of view by the atmosphere. Air- 
glow is emission from the recombination of electrons in 
the earth's atmosphere which were excited during the 
day by solar photons, and as such is a function of so- 
lar activi ty, time elapsed sin ce sunset, and geomagnetic 
latitude l)Leinert et alJ ll998). Airglow and atmospheric 
scattering vary throughout the night, moonlight varies 
from night to night, and zodiacal light varies from year 
to year. The background values from frame to frame cor- 
respondingly vary temporally by up to 10% throughout 
one run and 20% from year to year. 

Due to the temporal variations in the background, it 
is necessary to link the off-cluster backgrounds from ad- 
jacent frames to create one single background of zero 
counts for the entire cluster mosaic before averaging to- 
gether frames. To determine the background on each 
individual frame we measure average counts in approx- 
imately twenty 20 x 20 pixel regions across the frame. 
Regions are chosen individually by hand to be a repre- 
sentative sample of all areas of the frame that are more 
distant than O^/i^Mpc from the center of the cluster. 
This is well beyond the radius at which ICL components 
have been identified in ot her clusters ijFeldmeier et alJ 
120021: iGonzalez et al.l |2005). and is also beyond the ra- 
dius at which we detect any diffuse light in A3888. The 
average of these background regions for each frame is 
subtracted from the data, bringing every frame to a zero 
background. The accuracy of the background subtrac- 
tion will be discussed in |Q 

4.6. Extinction Correction 

After background subtraction, all flux in the frame 
originates above the atmosphere, and is subject to atmo- 
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spheric extinction (large angle scattering out of the line 
of sight). This is equally true of resolved sources and 
diffuse sources less than several degrees in extent. While 
extinction corrections are usually applied to individual 
resolved sources, that is not possible with the diffuse 
ICL. We correct entire cluster images for this extinction 
by multiplying each individual image by 10 7 "*/ 2 - 5 , where 
X is the airmass and r is the fitted extinction term from 
the photometric solution. This multiplicative correction 
is between 1.06 and 1.29 for the airmass range of our 
A3888 observations. 

4.7. Registration & Distortion 

To combine images, we align all 41 individual frames 
to one central reference frame. SExtractor positions of 
approximately 10 stars in each frame are used as input 
coordinates to the IRAF tasks GEOMAP and GEOTRAN, 
which find and apply x and y shifts and rotations between 
images. The GEOTRAN solution is accurate to 0.01 pixels 
(RMS). As an independent check of registration accu- 
racy, we confirm that the center coordinates of stars in 
the original images, as compared to the combined image, 
are the same to within 0.01 pixels. This uncertainty is 
negligible for our measurement which is made on much 
larger scales. In addition, the ellipticities of individual 
stars do not change with image combination, suggesting 
that no systematic errors in registration exist. Stellar 
ellipticities also show no variation across the frame, sug- 
gesting that there are no significant image distortions. 

4.8. Image Combination 

After pre-processing, background subtraction, extinc- 
tion correction, and registration, we combined the images 
using the IRAF routine imcombine with a rejection of 
— 3.5cr and +4.5<7. This range was chosen as a compro- 
mise between rejecting the cosmic rays (CRs) and allow- 
ing for some seeing variations in the peak flux of stars. In 
total, 16 and 25 900s exposures in the V— and r— bands, 
respectively, were averaged together. The final combined 
image is 4096 pixels (3.6hy C , Mpc) on a side. The central 
region (approximately l/i™ Mpc on a side) of the final 
combined V— band image is shown in Figure^ 

5. ANALYSIS 

5.1. Object Detection 

We use SExtractor both to find all objects in the com- 
bined frames and to determine their shape parameters. 
The detection threshold in both the V and r images was 
defined such that objects have a minimum of 6 contigu- 
ous pixels, each of which are greater than 1.5a above 
the background sky level. This corresponds to a mini- 
mum surface brightness of 26.0 mag arcsec -2 in V and 
26.4 mag arcsec -2 in r. The faintest object in the cat- 
alog has a total magnitude of 27.0 mag in V and 27.4 
mag in r, however we are complete only to 24.8 mag 
in V and 24.5 mag in r. We choose these parameters 
as a compromise between detecting faint objects in high 
signal-to-noise regions and rejecting noise fluctuations in 
low signal-to-noise regions. Shape parameters are deter- 
mined in SExtractor using only those pixels above the 
detection threshold. 

5.2. Object Removal & Masking 
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Fig. 1. — The central ~ IMpc (4.9 arcmin) of A3888 shown in 
the V band. The gray-scale is linear over the range 20.4-29.5 mag 
arcsec -2 . Note the three main groups of galaxies near the center 
of the cluster. A dozen objects in this image are stars; the rest are 
galaxies. 



To measure the ICL we remove all detected objects 
from the frame by either subtraction of an analytical 
profile or masking. Details of this process are described 
below. 

5.2.1. Stars 

Scattered light in the telescope and atmosphere pro- 
duce an extended point spread function (PSF) for all 
objects. To correct for this effect, we determine the tele- 
scope PSF using the profiles of a collection of stars from 
super-saturated 4th mag stars to unsaturated 14th mag- 
nitude stars. The radial profiles of these stars were fit 
together to form one PSF such that the extremely satu- 
rated star was used to create the profile at large radii and 
the unsaturated stars were used for the inner portion of 
the profile. This allows us to create an accurate PSF to 
a radius of 7', shown in Figure |2j 

The inner region of the PSF is well fit by a Moffat func- 
tion. The outer region is well fit by r -3 . There is a small 
additional halo of light at roughly 50 - 100"(200-400pix) 
around stars imaged on the CCD. Images of saturated 
stars located off of the field of view of the detector do 
not show this halo, indicating that it is due to a reflec- 
tion of light off of the CCD itself. We find that roughly 
1% of the total flux in the star is in this halo. There 
are 13 saturated stars within 3.8 Mpc of A3888 ranging 
from 11.6 - 15.2 V magnitudes. The nearest three satu- 
rated stars are 0.6, 0.8, and 1.0 /i^q Mpc from the cluster 
center and have 14.6, 13.4, and 11.6 V magnitudes, re- 
spectively. These stars do not directly contribute to the 
ICL measurement because they are not near enough to 
the center, do not have very bright magnitudes, and the 
PSF does not put very much power into the wings. We 
do a careful job of background subtraction, by tying to 
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Fig. 2.— The PSF of the 100-inch du Pont telescope at Las 
Campanas Observatory. The y-axis shows surface brightness scaled 
to correspond to the total flux of a zero magnitude star. The 
profile within 5" was measured from unsaturated stars and can be 
affected by seeing. The outer profile was measured from two stars 
with super-saturated cores imaged in two different positions on the 
CCD, on two different observing runs. The bump in the profile at 
100" is likely related to the position of the star in the focal plane. 
If the core of a star is imaged off the CCD, its profile does not show 
this feature, suggesting that the feature is caused by a reflection 
off the CCD itself. The outer surface brightness profile decreases 
as roughly r~ 3 , shown by the solid line. An r -2 profile is plotted 
to show the range in slopes. 

off-cluster flux, so that the PSF also does not affect the 
background measurement. 

For each individual, non-saturated star, we subtract a 
scaled r~ 3 profile from the frame in addition to masking 
the inner 30" of the profile (the region which follows a 
Moffat profile). Since we do not have accurate magni- 
tudes for the saturated stars in our own data, and to be 
as cautious as possible with the PSF wings, we have as- 
sumed the brightest possible magnitudes for these stars 
given the full USNO catalog errors. We then subtract a 
stellar profile with that magnitude and produce a large 
mask to cover the inner regions and any bleeding. We 
can afford to be liberal with our saturated star mask- 
ing since there are very few saturated stars and none of 
them are near the center of the cluster where we need 
the unmasked area to measure any possible ICL. 

5.2.2. Galaxies 

To make an ICL measurement we would ideally like 
to subtract a scaled analytical profile for each galaxy 
that would leave no residuals and would allow us to re- 
cover the area on the sky covered by cluster galaxies. We 
have attempted t o do this using 3 pu blicly av ailable algo- 
rithm s: GIM2D ijSimard et al.ll2002fl Galfit llPeng et all 
I2002|) . and the IRAF task ellipse l|Jedrzeiewskilll987j) . 
With these algorithms, we have employed a wide range 
of surface brightness profiles, including deVaucouleurs, 
Sersic, exponential profiles, and combinations thereof. In 




□ □H 

Fig. 3. — Images of observed galaxies, GALFIT models, and 
model-subtracted residual images are shown from left to right for 
3 different galaxies in A3888. All images are shown at the same sur- 
face brightness levels. The first row shows a fairly isolated galaxy 
in the outer regions of the cluster, which is well modeled by Gal- 
fit. The second and third rows show galaxies in increasingly denser 
environments, depicting well the limitations of galaxy modeling 
algorithms for galaxies in very dense regions. 



addition, we have used iterative techniques to alternately 
fit and remove galaxies in crowded fields. The technical 
challenges in fitting the galaxies, including galaxy de- 
blending, PSF effects and deconvolution, 2D profile fit- 
ting, and speed in performing many Fourier transforms 
have bee n previously discu ssed by several groups (see for 
cxa mple lPeng et al.ll2002l for a review ) . 

Figure shows representative results of modeling 3 
galaxies using Galfit: one isolated galaxy and two galax- 
ies in increasingly dense regions. These examples show 
that the algorithms perform well for isolated galaxies, but 
fail for galaxies near the core due either to difficulty in 
deblending many overlapping galaxy profiles or because 
the individual galaxies in such dense regions do not follow 
simple analytical profiles. It is not clear what the pro- 
files sho uld be of galaxies deep in the potential w ells of 
clusters llTruiillo et al.ll200HlFeIdmeier et al.ll2 004). The 
fact that A3888 is not a relaxed cluster clearly makes 
galaxy subtraction more difficult near the core than it 
would be in a CD cluster; A3888 has 3 main brightness 
peaks which contain 3, 7, and 12 galaxy cores in their 
densest regions, respectively. 

As it is not possible to cleanly fit the galaxies in this 
cluster such that the residuals (positive or negative) do 
not interfere with the ICL measurement, we have chosen 
to mask the galaxies. This gives us a well defined mea- 
surement of the ICL at the expense of forfeiting some 
area. Although we could model and subtract the more 
isolated galaxies in the outer regions of the cluster, it is 
in these regions that we can generously mask the galax- 
ies and still have enough pixels for an ICL measurement. 
Note that we do not replace masked pixels. Masked re- 
gions are simply removed entirely from the ICL measure- 
ment. 

We use the same masks for both bands so that all 
galaxies are masked to the same radius, thereby insur- 
ing a self-consistent measurement of the ICL color. We 
use the r— band image to define the masks as it has a 
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Fig. 4. — The fully masked, final V— band image of the central 
1.5 hyp 1 Mpc (7.3 arcmin) of A3888, smoothed to aid in visual 
identification of surfac e brigh tness levels. Masks are shown in their 
intermediate size, see H5.2.2I Large circular masks correspond to 
the locations of bright stars. The six gray-scale levels show surface 
brightness levels of up to 28.5, 27.7, 27.2, 26.7 mag arcsec -2 , and 
brighter than 26.2mag arcsec -2 . ELLIPSE isophotes are o verla id 
from 65 - 190". The tidal feature, C, shown also in Figure ITTI is 
clearly visible at center near the bottom of the image. 



deeper detection threshold (and thus larger detection ar- 
eas) than the V— band catalog. Objects are identified 
using SExtractor and masks are based on the isophotal 
detection area with a threshold of 26.4 mag arcsec _2 (1.5tr 
above sky). To be conservative in rejection, we scale the 
semi-major and semi-minor axes identified by SExtractor 
to increase the area of each galaxy mask by a multiplica- 
tive factor of 2 — 2.3, depending on the magnitude of the 
galaxy. To explore the effect of mask size on the profile 
shape of the ICL, we make two additional images with 
mask sizes that are 30% smaller and 30% larger than the 
original masks. We then measure the ICL three times 
with the three versions of mask sizes. Additional minor 
masking is done by hand to remove any remaining flux 
associated with resolved objects. These few regions are 
associated with small overlapping sources which are not 
correctly deblended by SExtractor. 

The total masked area within the central 1.2 Mpc of 
the cluster in each of the three mask sizes is 34%, 41%, 
and 49%. The masked fraction is much higher in the 
very center of the cluster and reaches nearly 100% in 
the inner 30 arcseconds. The increase in masked fraction 
is not directly proportional to the increase in mask size 
because the masks often overlap. Figure^lshows the final 
V— band image with intermediate-sized masks. 

5.3. Cluster Membership & Flux 

An interesting characteristic of the ICL lies in its com- 
parison to cluster properties including the cluster galax- 
ies themselves. We compare two methods below for mea- 
suring cluster membership and flux: (1) we identify mem- 
ber galaxies using our own 2-band photometry; and (2) 
we integrate the flux in a published galaxy luminosity 
function for this cluster. 
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Fig. 5. — The color magnitude diagram of galaxies in A3888. 
All galaxies detected in our data are plotted with gray symbols. 
Those galaxies which have membership information in the liter- 
ature are over-plotted with open triangles (members) or squares 
(non-members). The red sequence is clearly visible. Solid lines 
indicate a biweight fit to the red sequence with lc uncertainties. 



So me published veloci t ies are available in t he litera- 
ture i|Teague et al.lll99ft iPimbblet et al.ll2002fl and can 
be used to explicitly identify member galaxies. However 
these redshift surveys are not complete to our detection 
threshold, and can therefore not provide membership in- 
formation for all detected galaxies. Alternatively, we can 
estimate cluster membership using a color magnitude re- 
lation (Figure [SJ from our V and gunn-r images. There 
is a clear red sequence of galaxies where the brightest 
galaxies have V — r = 0.3±0.15. Those galaxies which 
lie within la of a biweight fit to the red sequence are 
taken to be cluste r members (functional form taken from 
iBeers et al.l ll990^ The slope of the red sequence is 0.1 
mags(color)/ 4mags(galaxy r magnitude). Those galax- 
ies which are redder than the red sequence are both gen- 
erally fainter implying that they are higher redshift back- 
ground galaxies and are not as concentrated toward the 
center of the cluster as all galaxies. The number of those 
very red galaxies per projected area is 38 ± 11% higher 
within 400kpc than without. Although some of these 
galaxies are undoubtedly members of the cluster, their 
spatial distribution does not allow us to make conclusive 
statements about their membership. Approximately 42% 
of the galaxies in the image are identified as members by 
this method. Of the galaxies with spectroscopically de- 
termined velocities, 78% of the 55 confirmed members are 
included in the cut; 54% of the 13 known non-members 
are also included. The red cluster sequence is a good tool 
for identifying clusters, but it is not a perfect method of 
determining membership as it is unable to cleanly distin- 
guish between member and non-member galaxies. 

We measure the total flux in all galaxies identified 
as members using corrected isophotal magnitudes from 
SExtractor. For these, SExtractor assumes a Gaussian 
profile to infer the flux beyond the isophotal detection 
threshold, corresponding to 26.0 V mag arcsec -2 and 
26.4 r mag arcsec -2 . As expected, the corrected mag- 
nitudes are brighter than the isophotal magnitudes by 
a full magnitude at the faint end of our detection limit. 
The total flux in galaxies within 700hy 1 kpc center of 
the cluster, as determined from the same galaxy catalog 
that was used for galaxy masking, is 3.9 x 10 12 L in 
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the V— band and 4.9 x 10 L Q in the r— band. We ex- 
pect the error on the total flux from this estimate to be 
greater than 30% which is mainly due to uncertainty in 
the membership determination. 

We can also determine cluster flux using the 
iDriver et al.1 l|1998j) luminosity distribution for this clus- 
ter, which is based on a statistical background subtrac- 
tion of non-cluster galaxies. It would be possible to do 
this with our own data, however lDriver et alJ i|1998f) have 
more uniform, large-area coverage to several magnitudes 
below M* at the redshift of the cluster. In addition, the 
authors pay careful attention to observing background 
fields which are up to 75' from the cluster center, at ap- 
proximately the same airmass, seeing, exptime, and UT 
as the cluster fields. Consequently, the background fields 
have the same noise characteristics and detection thresh- 
old as the cluster images, and sample the same large scale 
structures. They can therefore be used t o reliably deter- 
mine c ontamination of the cluster fields. IBernstein et alJ 
(1995) gives a careful account of the significant consider- 
ations in u sing this method, a ll of which are taken into 
account bv IDriver etafl l)1998|L 

We explo r e on e minor effect not discussed by 
IDriver et al.l Jj.998): the effect of gravitational lensing 
on the background galaxy counts. There are two com- 
peting effects which change the number and brightness of 
galaxies behind the cluster as compared to background 
galaxy counts in an off-cluster field. First, magnifica- 
tion of the background galaxies will artificially inflate 
the background counts behind the cluster, resulting in 
an underestimation of cluster galaxy flux. Second, all 
background objects behind the cluster will appear radi- 
ally more distant from the cluster center, which will ar- 
tificially decrease the background counts, resulting in an 
overestimate of the cluster galaxy flux. The change from 
an overall magnification to de-magnifi cation happens at 
z ~ 0.5. Following the method of iBroadhurst et all 
( 1995) to determine the strength of the de-magnification 
for A3888 at z = 0.15, we find a negligible degradation 
in the V— and i — band flux (< 0.2%), and therefore do 
not correct for it in the IDriver et al.l l|1998|) background 
c ounts. 

IDriver et alJ l(l"998j) use their i?-band luminosity dis- 
tribution to determine a dwarf to giant ratio, however 
we choose to fit it with a classical Schechter function 
(M* = -22.82 ± 0.28, a = -0.97 ± 0.09, xl = 0.71), 
which can then be used to determine a luminosity density 
for the cluster. We note that the luminosity distribution 
is not perfectly fit by a Schechter function at the bright 
end, due mainly to a small number of extremely bright 
galaxies, as is typical of clusters. Adopting a volume 
equal to that over which we are able to measure the ICL, 
1.4 Mpc 3 , and integrating the luminosity function down 
to very faint dwarf galaxies, = -11, the total lumi- 
nosity from galaxies in the cluster is 5.9 ±0.94 x 10 12 
in the R— band. Given galaxy colors from lFukueita et al.l 
( 1995), the total luminosity from galaxies in A3888 is 3.4 
±0.6 x 1O 12 L in V and 4.3 ±0.7 x 1O 12 L in the r-band. 
The difference between this value of total flux and that 
determined from our color-magnitude estimate of mem- 
bership is likely due to uncertainties in our membership 
identification and difference in detection thresholds of the 
two surveys. Although the two estimates are generally 
consistent, we adopt the total flux as derived from the 



luminosity distribution throughout the remainder of the 
paper. 

6. RESULTS 
6.1. Surface brightness profile 

After subtracting the stars and masking the galaxies, 
we fit the resulting image with the IRAF routine ellipse, 
a 2D, interactive, isophote fitting algorithm. Again, the 
masked pixels are completely excluded in this procedure. 
There are 3 free parameters in the isophote fitting: cen- 
ter, position angle, and ellipticity. We fix the center ( 
J2000.0, 22^34 m 26.0 s , -37°44' 07.2") and position angle 
(-70 degrees) to values found by ellipse based on the in- 
ner isophotes, and let the ellipticity vary as a function of 
radius. Fitted ellipticities range from 0.2 to 0.5. Allow- 
ing the center and position angle to vary results in worse 
fits. Stable fits are found from 60 — 250". From the fitted 
isophotes we identify a fairly smooth ICL profile over the 
range of 26 to approximately 29 mag arcsec -2 . The error 
on the mean within each elliptical isophote is negligible, 
as discussed in It is possible that the different seeing 
in the V and r— band images could unevenly affect the 
profiles. To address this issue, the V and r— band images 
have been convolved to the same seeing, and the surface 
brightness profiles re-measured. No significant change 
was found in the profiles. 

Note that we are not able to measure the ICL at radii 
smaller than 60" because that region is heavily masked. 
Most other ICL measurements focus on this inner region, 
leaving little overlap between this survey and previous 
work in other clusters. In clusters containing a cD galaxy, 
the diffuse component of the cluster has been found to 
blend smoothly into the cD envelope, and masking in the 
core of such clusters is not necessary (see most recently 
IGonzalez et al.ll2005lh 

We identify the surface brightness profile of the to- 
tal cluster light (ie., including resolved galaxies) for 
comparison with the ICL within the same radial ex- 
tent. To do this, we make a new "cluster" image, with 
color-determined, non-member galaxies masked out (see 
H5.3|) . A surface brightness profile of the cluster light is 
then measured from this image using the same elliptical 
isophotes as were used in the ICL profile measurement. 
This profile, in contrast to the ICL, is quite irregular, 
reflecting the clustering of galaxies. Substructure in the 
galaxy distribution is an indication of a young dynamical 
age for this cluster. 

Figure shows the surface brightness profiles of the 
ICL as well as the total cluster light as a function of 
semi-major axis in both the V— and j — bands. Results 
based on all three versions of mask size (as discussed in 
H5.2.2|l are shown. The uncertainty in the ICL surface 
brightness is dominated by the accuracy with which the 
background level can be identified, as discussed in J3 
Error bars in Figure show the cumulative uncertainties 
tabulated in Table 3. 

Two characteristics are evident from the surface bright- 
ness profiles. First, the inner region (200 — 400 h^kpc) 
has a notably steeper profile than the outer region. While 
the entire profile can be adequately described within the 
la uncertainties by a single exponential, a double expo- 
nential gives a better fit in the r— band (xt improves by 
50%) and a marginally better fit in the V— band. These 
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Fig. 6. — The surface brightness profile of the intracluster light 
as well as the total cluster light plotted as a function of distance 
along the semi-major axis in arcseconds. The axis at the top of 
the figure indicates corresponding physical scale in h^Mpc. We 
plot both V— and r— band data together for comparison. The 
bottom two lines on the plot are the ICL profiles; the r— band light 
is surrounded by solid shading and the V— band is surrounded by 
hatched shading. The shadings show the difference in ICL profiles 
produced by in creasin g or decreasing the area of the galaxy masks, 
as discussed in H5.2.2I The top two lines, without shading represent 
the total cluster light as measured in the same elliptical isophotes 
as the ICL; the dashed line represents the V— band light, the solid 
line represents the r— band. Also shown are the cumulative la 
errors for both bands as discussed in J7|and summarized in Table 
3. 




Fig. 7. — The V— band intracluster light and 2a error bars over- 
plotted with exponential fits. The best fit singlc(dashcd) and dou- 
ble (gray) exponential are shown. 



fits are shown in Figures [7| and [SJ We have also fit the 
ICL profile with DeVaucouleurs and Sersic profiles. Ac- 
ceptable fits can be found, however the best fit values 
are unphysical. Namely they have high exponents for the 
Sersic and unrealistically large effective radii for the De- 
Vaucouleurs profiles. The second general characteristic 
of the ICL is that it is more concentrated than the galax- 
ies, which is to say that the ICL falls off more rapidly 
with radius than the galaxy light. 

6.2. Spatial Distribution 

The ICL is aligned to within 10° of the position angle 
of the hot intracluster gas. Figure ED shows contours of 
XMM archival observations overlaid on our optical im- 
age. We interpret the alignment of the diffuse intraclus- 
ter light with the hot gas in the cluster as an indication 
that we are indeed measuring light which follows the 
gravitational potential of the cluster. In addition, the 
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Fig. 8. — Same as Figure l7l in the r— band. 




Fig. 9. — X-ray contours taken from XMM archival data are 
overlaid on our V— band optical image. Logarithmic contours are 
shown from 1 to 20 counts. The bright point source GOOh^ 1 kpc 
from the cluster center is a Seyfert I galaxy. 



ICL radial surface brightness profile is significantly dif- 
ferent than the galaxy surface brightness profile in both 
V and r, suggesting that the intracluster light component 
is at least in part distinct from the individual galaxies in 
the cluster. 

6.3. Color 

We measure an average V — r color of the ICL by bin- 
ning together three to four points from the ICL radial 
profile. Between 200 and 400hy 1 kpc, the innermost 



measured radii, the diffuse ICL has an average color of 
V - r ~ 0.3 ± 0.1. Beyond 400h^j L kpc the ICL becomes 
increasingly redder, such that by TOOh^Q 1 kpc the ICL has 
an average color oi V — r ~ 0.7 ± 0.4. The only char- 
acteristic color of the galaxies we have to compare with 
the ICL is the red sequence color (V — r = 0.3 ± 0.15). 
We have no definitive membership information for those 
galaxies off the red sequence. The color of the ICL in 
the inner 400h 70 1 kpc is roughly equivalent to the red 
ellipticals residing in the same part of the cluster, but 
significantly redder than several tidal features we detect 
(see The color of the ICL beyond 400hf 1 kpc is 

redder than the red sequence galaxies. The color of the 
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Fig. 10. — The ICL color versus radius in coarse radial bins based 
on the surface brightness in V and r as shown in Figure |S| Lower 
axis shows radius in arcseconds, upper axis shows radius in Mpc. 
The dashed line is the best fit linear function. The average color 
of the red cluster sequence and the tidal features are also shown 
for comparison 



diffuse ICL can be approximated as a simple linear func- 
tion of radius, with a slope of +0.1 per lOOhyQ 1 kpc and 
a y-intercept of —0.1. Figure ITUI shows the color profile 
and corresponding lex error bars. While this fit is clearly 
simplistic, the data do not warrant a more complicated 
fit. This red color gradient is opposite that which we ex- 
pect to find for the cluster galaxies. When looking at the 
color of galaxies as a function of distance from the center 
of the cluster, we find a flat or slightly blueward profile 
such that the galaxies get slightly bluer with increasing 
radius. Therefore the ICL color profile is distinct from 
the galaxy color profile. 

Using the p op ulatio n synthesis models of 
IBruzual fc Chariot! |2003) we can obtain rough 
constraints on the age and metallicity of the stel- 
lar population which contributes to the ICL. Because 
the total range in color is not large and because the 
age-metallicity degeneracy limits our conclusions, we 
limit our discussion to two regions, r < 400hf 1 kpc and 
r > 400h 70 1 kpc, rather than individual points along the 
full radial profile. The stellar evolution models begin 
with a starburst of user defined strength and extent, a 
Salpeter initial mass function, and a standard ACDM 
Universe. The stars then evolve along the Padova 1994 
spectral evolution tracks. Within this scheme, the 
simplest scenario is an instantaneous starburst with a 
single formation epoch and metallicity. For this case, 
the red color of the ICL in the outer 400 — 700/i 70 1 kpc, 
V — r ~ 0.6, is consistent with a stellar population which 
formed at redshifts 1 <z/ < 10 (7 - 13 Gyr ago) with an 
initial abundance of 1.0Z© < Zicl < 2.5Z©. The color 
of the ICL in the inner 200-400hf 1 kpc, V - r ~ 0.3, 
allows the minimum age of that range to be lowered, 
where the most recent allowable formation is ~ 5 Gyr 
ago (zf < 1) with an initial abundance of 0.2 - 0.5 solar. 
Allowing an extended burst of duration 10 to 100 Myr 
has a minimal effect on the V — r color. Allowing an 
exponentially decaying star burst with an e-folding time 
of lGyr, the population becomes overall 0.02 - 0.06 
magnitudes bluer, depending on the initial formation 
redshift. For the ICL in A3888, an exponential star 
formation history therefore implies even higher metal- 



licities or earlier formation. Finally, simulations with a 
constant star formation rate of 1M Q /yr create very blue 
stars. It is not possible to form a stellar population with 
a constant star formation rate which has V — r = 0.6. 
Implications of these models are discussed further in 

El 

6.4. Fractional Flux 

The ratio of ICL flux to total cluster flux can help 
constrain the importance of galaxy disruption in the 
evolution of clusters. To identify the total flux in the 
ICL, we integrate the single exponential fit to the ICL 
surface brightness profile (see i)6.1f) over the range 0- 
200" (~ TOOhyp 1 kpc). As we are not able to measure 
the ICL at radii smaller than ~ 60", this requires an ex- 
trapolation into the center of the cluster. Note that the 
single exponential fit, which is dominated by the slope 
of the ICL profile at larger radii, gives significantly less 
light in the core than the double exponential fit (see Fig- 
ures Hand |HJ, and is therefore a conservative estimate of 
ICL flux. The total flux in the intracluster light is then 
4.5 ± 1.3 x 10 11 Lq in V and 5.9 ± 2.2 x 10 n L Q in r, 
where these errors are the full errors as described in SJ3 
This value is equivalent to the full disruption of roughly 
7L* galaxies. 

We consider 4 modifications to this estimate of the to- 
tal ICL flux. First, we consider a correction for that 
volume of the cluster which is filled with galaxies, since 
no ICL can exist in that volume. While lines of sight 
intersect galaxies over most of the area near the center 
of the cluster, the galaxy filling factor is less than 3% by 
volume, even inside 200hy 1 kpc (60 arcsec, projected). 
So it is reasonable to assume that intracluster stars do 
exist in that volume and we need make no correction for 
the filled volume. Second, we can determine a hard lower 
limit to the ICL flux by assuming that there is no ICL in 
the inner 60". This correction, albeit extreme, would de- 
crease our estimate of the total ICL flux by 30%. Third, 
we make a less extreme correction by assuming a flat 
core region instead of th e exponential extrap olation. A 
flat profile is suggested bv lAguerri et alJ l)2005t) , although 
those results are in Virgo, where the center of the ICL 
is not defined and the measurement is based on small 
area coverage, which does not allow for an elliptical pro- 
file determination. A flat core region would decrease our 
estimate of the total ICL flux by 5%. Fourth, we con- 
sider low surface brightness galaxies below our detection 
threshold which could contribute to, and therefore be an 
error in, the inferred ICL flux. To account for these very 
faint galaxies, we integrate the cluster galaxy luminosity 
function from our detection limit (Mb = -15.22) to Mr 
= -11.0. Due to the extremely low detection threshold of 
this survey (7.6 magnitudes below M*), and the appar- 
ently flat faint end of the luminosity function (a = —0.97, 
see H5.3(l . only 0.07% of the total galaxy flux could come 
from galaxies this faint. As this contribution is not sig- 
nificant, we make no correction for this effect. 

Adopting the total galaxy flux found from the lumi- 
nosity function in H5.3I we find that the ICL accounts for 
13 ± 4% of the total U-band cluster light and 13 ± 5% of 
the r— band cluster light within 700h 70 kpc of the center 
of the cluster. The range in these values comes from the 
combination of all uncertainties in the ICL measurement 
coupled with the uncertainty in the total cluster flux, as 
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Fig. 11. — Three low surface brightness features from final 
V— band image. The top two panels show features located near 
the center (features A & B in Table 2). The bottom image shows 
the larger feature (C) located SOOh^, 1 kpc from the center of the 
cluster. 



discussed in SJ3 The galaxy light and the ICL decrease 
with radius, but since we don't accurately know the slope 
of either of them at large radii, we compare fluxes within 
the same volume over which we have reliable data. This 
fraction is only relevant at this radius, and is likely to be 
lower when taking into account the entire virial radius of 
the cluster, since the ICL is centrally concentrated and 
not evenly distributed throughout the cluster. On the 
same note, ICL measurements at smaller radii are likely 
to find a higher fraction of the total flux in the ICL be- 
cause of the steep ICL profile, and because the volume 
involved is much smaller. For example, if we assume we 
can only measure the ICL in the inner 600hy 1 instead of 
700h 70 1 kpc, we find a fractional flux of 19% in both V 
and r, an almost 50% increase over the measured 13%. 

6.5. ICL Substructure 

Using the technique of unsharp masking, we find 3 pos- 
sible tidal features, all within the central SOPh^ 1 kpc of 
A3888, identified as A, B, and C in Figure [n] Arcs 
A and B are both roughly lSh^kpc xSh^kpc and are 
near to the center of the cluster. Arc C is a diffuse, 
tail- like feature at SOOh^kpc from the center and covers 
lSOh^kpc x20h 70 1 kpc (see Table 2). All three features 
are blue, V — r ~ 0.0, with a combined flux equivalent to 
one r — 20.8 magnitude galaxy (0.1M*). These objects 
are unlikely to be gravitational arcs since they are not 
oriented tangentially to the cluster potential. 

The diffuse nature of the large feature, C, sug- 
gests that it is tidal. This object is very sim- 
ilar to the large arcs fo und in Coma, Centaurus, 
and A1914 iGregg fc Westjll998t ICalcaneo-Roldan et all 
120001 iFeldmeier et al.ll2004j) "which are included in Ta- 
ble 2 for reference. In general it is of similar size to, 
but slightly fainter than, those found in the other clus- 
ters. It has slight curvature and appears to connect 



to a pair of galaxies (left side of the image as dis- 
playe d) that could be in the mid st of an interaction. 
Both l Calcaneo-Roldan et alJ l)2000|) and IFeldmeier et alJ 
(2004) find through numerical simulations that these 
types of arcs are typical of recent tidal interactions be- 
tween luminous spiral galaxies and massive cluster ellip- 
ticals. Spectroscopy to confirm its origin at this faint 
surface brightness is not currently possible. 

We cannot rule out the possibility that the 2 smaller 
features in our 0.06 degree -2 of cluster imaging are low 
surface brightness (LSB) galaxies seen edge-on. In field 
surveys, surface densities of the dimmest LSB galaxies 
(23 < no < 25V mag arcsec -2 ) a re at least 0.01 galaxies 
degree -2 (|Dalcanton et alJ ll.997). In clusters, although 
there are overall higher space densities of galaxies, LSB 
galaxies run the risk of getting tidally disrupted in the 
harsh environment of cluste r centers. In a surve y of the 
Cancer and Pegasus clusters. ICTNeil et all ltl997h find 1.6 
galaxies per square degree with central surface brightness 
dimmer than 21.2 in V. The 2 candidates with aver- 
age surface brightnesses of roughly 25mag arcsec -2 in V 
in this cluster represent a higher density than found in 
these surveys. In addition, they do not have clear cen- 
ters. Both of these facts suggest that they are not LSB 
galaxies. However, it is likely that even if these are LSB 
galaxies, they will not remain bound systems for long 
in the high density cluster environment, and we there- 
fore consider them to be contributors to the ICL in the 
following calculation. 

We briefly examine the importance of all three tidal 
features in contributing to the ICL over a Hubble time 
to see if they can account, in whole or in part, for the 
ICL found in the cluster. Cluster crossing time is esti- 
mated to be 4.5 Gyr given a virial radius of 3.7 Mpc an d 
a temperature of lOKeV ijReiprich fc Bohringerl 12002). 
We assume both a constant rate of formation, and dissi- 
pation of tidal features in approximately 1 crossing time. 
From this we conclude that in 1 Hubble time, approxi- 
mately one half of an M* galaxy will be contributed to 
the ICL through the visible tidal features such as these. 
This simple calculation suggests that these features can- 
not account for the current ICL flux, however it is feasible 
that there was a variable interaction rate in the history 
of this cluster. Further substructure could also be hidden 
below our surface brightness detection threshold. 

At the distance of A3888, the flux of a single glob- 
ular cluster (My > 10 mag) spread over one seeing 
disk (3.5/iY 1 kpc)is many magnitudes below our sur- 
face brightness detection threshold. Therefore we are 
not sensitive to intracluster globular clusters which 
have been studied b y other groups in n e arby c l usters 
llordan et alJ 120031 Bassino et all 12003 iHilkerl 1200.1 
Marin- Fra nch Sz, Apariciol l2003) . 

6.6. Group 

In addition to the main cluster ICL, we detect ex- 
cess diffuse light around a group of galaxies which 
are 1.35hf x Mpc from the center of A3888, (J2000) 
22 h 34 m 48.5 s , -37°39' 19.58". There are two galaxies 
centered in this diffuse component, separated by only 2". 
The spatial extent of the group appears to be 200h 70 1 
kpc, within which there are 60 galaxy peaks detected 
by SExtractor. Independent of the ICL component, the 
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TABLE 2 

Substructure in A3888 



Object 


radial distance 
h^kpc 


length 
hf^kpc 


width 
h^kpc 


V -r 


M v 




mag arcsec 2 


H(r) 
mag arcsec - 2 


A 


66 


17 


5 


-0.05 


-18.5 


-18.5 


24.6 


24.5 


13 


155 


15 


5 





-18.1 


-18.1 


25.3 


24.9 


C 


720 


132 


20 


0.03 


-20.5 


-20.5 


26.4 


26.7 



Coma 1 100 130 15-30 0.57(R) -18.8 -19.4(R) 26.9 26.3(R) 

Centaurus 2 170 171 1 ■ ■ • ■ • • -14.8(R) ■ ■ • 26.1(R) 

A1914 3 75 160 30 ■■• ■•• ••■ 26.1 



Note. — 1: IQregg fc West! I199SD 2: ICalcaneo-Roldan et all 1200(11 3: IFeldmeier et alJ i2004fl 



group is identified in the density distribution of cluster. 
Velocities are available only for the central galaxies in 
the group, however these suggest that the group is co- 
spati al with A3888 ijPimbblet et al.ll2002t iTeaeue et all 
1990). Within the 200hf 1 kpc extent of the group, we 
find 1.7 x 10 10 L Q in V and 2.6 x 10 10 L Q in r above 
background, which is equivalent to approximately 20% 
of an M* galaxy. 

The average color of this diffuse component based on 
total flux within 200hy 1 kpc is V — r = 0.5, which is 
again redder than the cluster galaxies and consistent with 
the color of the main cluster ICL at large radii. The 
accuracy of the fluxes and hence the colors is limited by 
the accuracy in masking since it is a simple sum over 
the pixels in the group region. We estimate the error 
in masking to be less than 30% based on our work with 
varying the mask size (see i)6.1l and Table 3). 

This secondary ICL concentration is consistent with 
galaxy interactions and ram pressure stripping occur- 
ring in an in-falling group ("pre-processing"). Such 
pre-processing has been shown in simulations to affect 
gala xies before they fall into the main cluster poten- 
tial l|Willman et all 120041 iFuiital 12004 and references 
therein). This is also consistent with recent measure- 
ments of a small amount of ICL in isolated galaxy groups 
lICa.stro-Rodrfguez et aljEfjoll IDurrell et alJEfiol. 

7. ACCURACY LIMITS 

The accuracy of the ICL surface brightness is limited 
on small scales (< 10") by photon noise. On larger scales 
(> 10"), structure in the background level (be it physical 
on the sky or instrumental) will dominate the error bud- 
get. We determine the stability of the background level 
in the image on large scales by first median smoothing 
the masked image by 75". We then measure the mean 
flux in thousands of random 1" regions more distant than 
0.8 Mpc from the center of the cluster. The standard de- 
viation of these regions is 29.5 mag arcsec -2 in V (0.06% 
of sky), and 28.8 mag arcsec -2 in r (0.01% of sky). His- 
tograms with gaussians overlaid are shown in Figures IT2l 
& EH The histograms are not perfect gaussians. This 
is likely due to the fact that the background level in- 
cludes both a symmetric gaussian and positive sources 
which are below the detection threshold. The offset of the 
gaussian portion of the histogram represents the statis- 
tical difficulty in measuring the mean value of the back- 
ground in any one image. Regions from all around the 
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Fig. 12. — Histogram of the mean values in 1" X 1" regions in 
counts in the fully masked V— band image. All regions are greater 
than 800h^Q kpc from the center of the cluster. A gaussian fit to 
the data is overlaid. 
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Fig. 13. — Same as Figure [T2l in the r— band. 



frame are used to check that our accuracy limit is uni- 
versal across the image and not affected by location in 
the frame. This empirical measurement of the large-scale 
fluctuations across the image is dominated by the instru- 
mental flat-fielding accuracy, but includes contributions 
from the bias and dark subtraction, physical variations in 
the sky level, and the statistical uncertainties mentioned 
above. 

This empirical measurement of the large-scale back- 
ground fluctuations is likely to be a worst-case esti- 
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TABLE 3 
Error Budget 



Source 



la uncertainty 
(V) (r) 



contribution to ICL uncertainty (%) 
jU(0"- 100") n(!00"- 200") total ICL flux 
(V) (r) (V) (r) (V) (r) 



background level" 


29.5 mag arcsec 2 28.8 mag arcsec 2 


14 


18 


39 


45 


24 


31 


photometry 


0.02 mag 0.03 mag 


2 


3 


2 


3 


2 


3 


masking 


variation in mask area ±30 


5 


5 


14 


19 


9 


12 


std. dev. in mean c 


32.7 mag arcsec -2 32.7 mag arcsec -2 


3 


2 


2 


1 


3 


1 


(total) 




15 


19 


41 


50 


26 


33 



cluster flux <i 



16% 



16% 



Note. — a: Large scale fluctuations in background level arc measured empirically and include instrumental calibra- 
tion uncertainties as well as and true variations in background level (see [J^J- b: Object masks were scaled by ±30% in 
area to test the impact on ICL measurement (see H5.2.2I . c: The statistical uncertainty in the mean surface brightness 
of the ICL in each isophotc. d: Errors on the total cluster flux arc based on errors in the fit to the luminosity function 
(sec JfT3l . 



mate of the accuracy with which we can measure surface 
brightness on large scales because it is derived from the 
outer regions of the image where only 6-10 individual 
exposures have been combined. In the central regions 
of our imaging (r < BOOh-Tp 1 kpc) , roughly twice as many 
dithered images have been combined, which has the effect 
of smoothing out large-scale fluctuations in the illumina- 
tion pattern to a greater degree. We therefore expect 
greater accuracy in the center of the image where the 
ICL is being measured. 

All sources of uncertainty are listed in Table 3. In 
addition to dominant uncertainty from our ability to 
measure the large-scale fluctuations on the background 
as discussed above, we quantify the contributions from 
the photometry, masking, and the accuracy with which 
we can measure the mean in the individual elliptical 
isophotes. In total the error on the ICL flux is 26% in 
V and 33% in r, which in addition to a 16% error in the 
total cluster flux, leads to a 30-40% uncertainty in the 
fractional flux. 



8. DISCUSSION 

We measure a diffuse intracluster component in A3888 
to a radius of ~ TOOhyg 1 kpc in the V— and r— band down 
to 28.9 mag arcsec -2 and 28.2 mag arcsec -2 respectively. 
We discuss here the physical implications of color, total 
flux, and profile shape of the ICL. 

8.1. Color Implications 

Color information can place constraints on the age 
and metallicity of the progenitor population of the ICL, 
thereby shedding light on the dominant physical mech- 
anisms and timescales for galaxy disruption. Color in- 
formation may also be able to differentiate between the 
morphological types of progenitor galaxies. 

The color of the ICL in A3888 is consistent with some 
previous observational res ults in othe r clust ers, although 
those results vary widely. iSchomberU l)1988l) and iMackiel 
(1992) have found a wide range of results for cD en- 
velopes, from blue to red, with and without color gradi- 
ents. These surveys have typically been sensitive to a dif- 
fuse component within lOO-lSOhyQ 1 kpc, a much smaller 
radial extent than this survey. Recently, iGonzalez et alJ 



(2000) have found a mild color gradient where the ICL 
becomes redder with radius by W — I = 0.25 ± 0.08 from 
10 to 70hjQkpc. A gain, our observations cover a much 
larger radial region of A3888 and we have no information 
on the ICL in the core region because it contains several 
complicated, unmerged clu mps. Over a s i milar range in 
radius to our measurement, Zi betti et alJ l)2005[K from a 
stack of hundreds of SDSS clusters, find an ICL includ- 
ing the BCG which is similar in color to the galaxy light 
and has a flat or slightly blue color gradient with radius. 

If the ICL is composed of stars stripped from galaxies, 
its color relative to the galaxies is indicative of the epoch 
when it was stripped in the following sense. If the ICL 
is redder than the cluster galaxies, it is likely to have 
been stripped from the galaxies at early times (higher 
z). Stripped stars will passively evolve toward red col- 
ors, while the galaxies will continue to form stars. If, on 
the other hand, the ICL is of similar color to the average 
galaxy, the ICL is likely to have formed fr om the ongo- 
ing st ripping of stars (via harassment as in iMoore et alJ 
fl996^ . In this case the stripped stars should have roughly 
the same color at the current epoch as the galaxies at the 
current epoch. This picture is complicated by the fact 
that clusters are not made up of galaxies which were all 
formed at a single epoch and that we don't know the 
star formation rates of galaxies once they enter a cluster. 
While these simple trends hold for the colors of intra- 
cluster stars compared to galaxies, the color difference 
between passively evolving stars and low star forming 
galaxies may not be large enough to detect. 

Cluster evolution is complex due to a myriad of envi- 
ronmental influences. Several groups have produced hi- 
erarchical, ACDM simulations of clusters which include 
radiative cooling, star formation, and various feedback 
mechanisms, but differ primarily in star formation pre- 
scriptions and numerical resolution. These models can 
be divided based on their broad, empirical predictions 
for the color /for mat ion epoch of the ICL. 

Theoretical models in which the ICL forms early in the 
cluster history all suggest an ICL which is older, redder 
than the galaxy population because the galaxies continue 
to form new stars and the refore have younger mean ages 
than the ICL population (Dubinski 1998;: lMurante et alJ 
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ommer-Larsen et al. 2005). This is generally con- 
sistent with our results in the outer re gions of A3888. 
Specifically. ISommer-Larsen et alJ l)2005f) predict a slight 
color gradient in B — R such that the ICL becomes 0.1 
magnitudes bluer from to BOOh^, 1 kpc, while our data 
suggest the opposite trend with radius in V — r. 

Theoretical models in which the ICL forms throughout 
the cluster lifetime generally pre dict a younger , bluer 
intracluster stellar population llWillman et al] 120041 
IMoore et a.l.l[l999t IGnedinll2003T: IBekki et alJ !2001) since 
more recent stripping will have the chance to pull newly 
formed stars out of galaxies. Ongoing stripping is con- 
sistent with our results within about 400h^ 1 kpc, where 
the ICL is roughly the same color or slightly bluer than 
the r ed cluster sequ e nce. In an N-body + SPH simula- 
tion, iWillman et~afl l|2004D find that 50% of intracluster 
stars come from M* or brighter galaxies, which means 
the color of the intracluster stars should be in accord 
with the color of the outskirts of bright cluster galax- 
ies or equivalently the color of intermediate luminosity 
galaxies. The intermediate luminosity galaxies which we 
consider members of A3888 have a color in the range of 

V — r = 0.1 — 0.4. Our results are consistent with this 
prediction in the inner regions of the cluster. 

Recent observations of some int racluster HII regions 
(such as those found in Virgo bv IGerhard et all 120021: 
IRvan- Weber et all 120041: ICortese et all I2004D indicate 
that it may be possible for some intracluster stars to 
form in situ. In this case, the ICL color will still depend 
on the formation epoch. If the ICSP is just now form- 
ing in int racluster HII regions (such as those found in 
Virgo bv IGerhard et aLll2002t IRvan- Weber et al.ll2004t 
ICortese et alJ I2004T) then it will be blue, however if it 
formed earlier in cluster formation, then the ICSP will 
passively evolve toward redder colors. Only if the ICL 
were significantly bluer than the existing cluster galaxies 
could it be possible to definitively state from the color 
that a significant fraction of the ICL formed in situ. Since 
this is not the case in A3888, our results cannot constrain 
the formation site of the ICL. 

U sing N-body s imulations, both IMoore et alJ ([1999) 
and iGnedinl ((2003) find that low density galaxies (LSB 
and dwarf galaxies) are the main contributers to the ICL. 
LSB galaxies in Sloan data have a color range of B — 

V = 0.4-1.3 ijKniazev et al.ll2004|) which corresponds to 

V — r = —0.2 — 0.7 . Dwarf g alaxies in Coma have V — 
r ~ 0.5 ±0.3 (|Trcntham 1998). This range is sufficiently 
broad that it is consistent with the ICL at all radii in 
A3888, implying that the ICL could have origins in LSB 
or dwarf galaxies. 

In summary, the ICL in the outer regions of A3888 is 
consistent with the predictions for a stellar population 
which formed at redshifts higher than 1 and is signif- 
icantly metal-rich, implying an ICL which forms early 
with the collapse of the main cluster. The ICL in the 
center of A3888 (r< 400h7 1 kpc) is consistent with pre- 
dictions for a relatively younger population. This implies 
that within some core radius harassment type interac- 
tions are the dominant mechanism. 

8.2. Fractional Flux Implications 

Another clue to the dominant mechanism driving evo- 
lution in clusters comes from correlating ICL properties 
with the properties of the parent cluster. For example, 



a trend in ICL fraction with cluster mass but not red- 
shift, richness, or morphology would indicate that mass 
was the dominant mechanism which could predict ICL 
fraction. The calculation of the fractional ICL flux de- 
pends on many observational parameters including the 
surface brightness and radial limit of the ICL measure- 
ment itself; the surface brightness at which individually 
bound, resolved sources are distinguished from the ICL; 
and the volume over which the ICL flux and galaxy flux 
are measured. As these parameters vary widely in work 
by previous groups, it is difficult to make meaningful 
comparisons with results for other clusters in the litera- 
ture. In addition, A3888 is a very massive cluster, and 
is not as simple as clusters with a cD or clear bright- 
est cluster galaxy (BC G). A1914 is the only c luster with 
an ICL measurement (|Feldmeier et alJ 1200 4') which has 
overall similar characteristics to A3888. With similar de- 
tection limits to those employed here, those authors find 
an ICL fraction of 7% in the V— band, which is generally 
consistent with our results for A3888. 

With these observational complications and cluster 
parameters in mind, we can only generally conclude 
that previous measurements of the ICL in clusters over 
a wide range in redshift (0.003 < z < 0.41) and 
mass (1 - 35 x 10 14 M Q ) are roughly 10%. There 
are no obvious trends with mass or redshift, although 
there are some notew orthy outliers at 50% for Coma 
llBernstein et al.l fl995h and 0% for A1689 <)Gudehusl 
I1989|) . The IBernstein et akl ljl995|) result covers a small 
radial extent and is therefore biased toward higher 
fractiona l flux values. It is difficult to interpret the 
iGudeh us (1989;) measurement due to disparate meth- 
ods. There is a some evidence that the ICL fraction 
is depen dent upon cluster mor phology; B/M type I 
clusters JTheuns fc Warredll997t iFeldmeier et al.ll2002l: 
lUson et "aTTIT99*1iriVikh'e^ have a re- 

ported average ICL fraction which is marginally higher 
than t hose with B/M type II I llVilchez- Gomez et alJ 
1991 iFeldmeier et al] l200l lArnabqldi et alJ 120031: 



Ferguson et aljl!998l IDurrell et al.ll2002^ ~ However poor 



morphological classification, small number statistics, and 
widely disparate methods and accuracies among the dif- 
ferent measurements make any possible trends difficult 
to quantify. 

In comparing the observed ICL with simulations, it is 
important to note that the simulations generally report 
the fractional light in the ICL out to much larger radii 
(e.g., r2oo) than its surface brightness can be measured 
observationally. At smaller radii, the predicted ratio of 
ICL t o galaxy light would be larger. Bearing this in 
mind, Willma n et al] l|2004j) finds a lower limit for frac- 
tional flux in the ICL of 10-22% at r2oo for a Virgo-like 
cluster from z = 1.1 — (increasing fractional flux with 
time). At the maximum radius of our ICL measurement 
(~ 0.3r2oo), the fraction would presumably be higher by 
at least a factor of two, making it larger than we observe 
in A3888. Other predictions are similarly high. For a 
cluster with the ma ss o f A3888 (25 x 10 14 Kn M ), both 
iLin fc Mohrl l)2004|) and lMurante et al] l)2004|) predict an 
ICL fraction in excess of 40%. To be consistent with their 
predictions, this cluster would require a factor of greater 
than 100 lower mass to have only 10% ICL, and although 
this cluster is not dynamically relaxed, such large erros 
in mass are not realistic. 
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In summary, we find an ICL fraction which is roughly 
compatible with observed ICL measurements in other 
clusters. However, our measurement differs signifi- 
cantly from theoretical estimates, particularly consider- 
ing A3888's large total mass. The dynamical state of 
A3888 may contribute to this discrepancy, as morphol- 
ogy may have a significant influence over ICL fraction. 
We emphasize again that A3888 is not a relaxed cluster; 
it does not have a cD galaxy and it's X-ray isophotcs 
are not circular. If we have caught this cluster as it 
is just now entering its major merger phase, we would 
expect a low ICL fraction as compared with a cluster 
at this redshift, mass, and richness which had already 
reached dynamic equilibrium. In contrast, when examin- 
ing Coma, a cluster with an extremely high ICL fraction 
but lower mass, we note that it's morphology indicates 
that it has already undergone significant merging to pro- 
duce 2 cD-like galaxies. If morphology is the dominant 
influence on ICL flux, we should find A3888 to have a 
similar ICL fraction to other clusters with similar mor- 
phologies. Comparable measurements of the ICL in our 
remaining sample will help to resolve this issue. 

8.3. Profile Shape Implications 

As discussed in tJ6.ll the profile of the ICL is generally 
stepper at smaller radii. In particular, the steepening 
profile near the core region of the cluster is associated 
with the 3 apparently merging groups of galaxies in the 
center of the cluster. The recent interactions in the cen- 
ter have likely added and continue to add ICL, which 
is likely to eventually relax into a BCG and BCG halo. 
The profile in the outer region of the cluster is consistent 
with previous measur ements of BCG envelo pes which fol- 
low shallower pr o files I I Gonzalez e \ aL^QQ a). In a ddition, 
iBernstein et~aLl (|1995h and iZibetti et alTl|2005H find a 
steeper profile for the Coma cluster and for a stacked 
profile of hundreds of Sloan clusters. This steeper profile 
is consi stent with some theor etical predictions, particu- 
larly bv lMurante et all (|2004f) based on a hydrodynami- 
cal simulation including radiative cooling, star formation, 
and supernova feedback. 

9. CONCLUSION 

We have presented results for the first of ten clusters 
in our sample. We have identified an intracluster com- 
ponent in A3888 to ~ TOOh^Q 1 kpc from the center of the 
cluster down to 28.9 mag arcsec -2 in the V— and 28.2 
mag arcsec -2 in gunn-r band. This ICL component is 
aligned with the hot gas in the cluster, which is evidence 
of its correlation with the underlying mass distribution. 
There is a second diffuse component around a group of 
galaxies l.Lhy^Mpc from the center of the cluster which 
is consistent with pre-processing in an in-falling group. 
In addition to these two diffuse ICL components, we find 
3 low surface brightness features consistent with being 
remnants from tidal interactions. 

Beyond 400h 70 1 kpc from the center of the cluster, the 
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ICL is redder than the galaxies, implying an older popu- 
lation of stars. Inside of 400h 70 1 kpc the ICL has a similar 
color to the galaxies. We interpret this color gradient in 
the ICL (y — r = 0.3 — 0.7, from inner to outer) as 
evidence of younger intracluster stars in the center of 
the cluster. Consequently, we suggest that more than 
one process is likely stripping stars from cluster galaxies. 
Specifically, harassment type interactions are still ongo- 
ing in the center of the cluster while galaxy mergers may 
have played a significant role earlier in the history of the 
cluster. 

We find that the ICL component in A3888 does not fol- 
low the same light profile as the resolved sources, but has 
a smoother and slightly steeper profile than the galaxies. 
Due to a steepening profile within 400h^ 1 kpc, the ICL 
profile can be described by a double exponential func- 
tion. A double profile is consistent with ongoing dynam- 
ical activity in the center of this cluster producing a new 
population of intracluster stars. 

Comparing the ICL to cluster galaxy flux, we find that 
the ICL component in A3888 accounts for roughly ql3% 
of the total cluster flux within TOOh^ 1 kpc (~ 0.3r 2 oo)- 
This value is low compared to the theoretical predictions 
for a cluster of this mass, and may be partly due to the 
fact that A3888 appears to be a dynamically young clus- 
ter. The ICL in A3888 will likely increase with time due 
both to contributions from an in-falling group as well 
as through major mergers in the center to create a cD 
galaxy. 
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Table 1: Cluster characteristics 



cluster ra (J2000) dec (J2000) 1 b 

name (hr:mm) (dd:mm) (degrees) (degrees) 



B/M Richness L x (0.1-2.4 keV) 
Class 10 44 ergs/s 



(7„ Number of n(HI) 1 

km/s confirmed members 10 20 cm~ 2 



lensing 
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A4010t 
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A3984t 
A014lt 
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23:57 
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23:13 
23:31 
22:34 
23:15 
01:06 
00:14 
22:58 



-34:40 
-30:34 
-28:52 
-21:38 
-36:30 
-37:43 
-37:48 
-24:35 
-30:2 
-34:47 



356.8 
17.99 
19.46 
41.37 
359.0 
3.96 
359.0 
175.3 
8.90 
8.32 



-76.06 0.048 r 
-58.50 0.058 r 



I 

II 



-80.98 0.062 r III 

-66.97 0.087 r II-III 

-70.60 0.096 r I II 

-59.40 0.15F III 

-67.19 0.18F II-III 

-85.93 0.23 r III 

-81.24 0.308 r III 

-64.78 0.31 a II-III 



3.09' 
1.86 1 
2.55 1 
2.47 1 
5.55 1 
14. 52 1 
9.18 1 
12. 62 1 
22.05 1 
7.0 U 



8451?™' 

g27 +120 m 

628±j£ m 
12471'"' 

no2t;^° 
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-2021, 
201 



1388±^ 8n 



4 5 d.p 

122 q - f ' k 
127i.P.k 

5 e,o.b 

30>- p 

69 s 

f 

a 

363 c - s ' h 
380 gJ ' h 



1.1 
1.1 

1.8 
2.0 
1.4 
1.2 
1.7 
1.8 
1.7 
2.0 



Dahle et al. (2002) 
Smail et al. (1991) 
Smail et al. (1991) 



O 
t" 1 



Notes: f: clusters for which we have done a photometric and spectroscopic survey for additional membership information (see §2) a: Abell, Corwin, k Olowin (1989). b: Batuski et al. (1999). c: 
Busarello et al. (2002) d: Chen et al (1998). e: Ciardullo, Ford, k Harms, (1985). f: Collins et al. (1995). g: Couch k Newell (1984). h: Couch k Sharpless (1987). i: Couch et al (2001). j: den 
Ilartog (1995). k: De Proris et al. (2002). 1: Ebeling et al. (1996). m: Girardi et al. (1998). n: Girardi k Mezetti (2001). o: Kowalski, Dlmer, k Cruddace, (1983). p: Mazure et al. (1996). q: 
Stein (1996). r: Struble k Rood (1999). s: Teague, Carter, k Gray, (1990). t: Wu, Xue, k Fang (1999). u: De Filippis et al. (2004). 



